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We study the upper mass bound of the lightest neutral Higgs scalar in the NMSSM 
using the RGE analysis. We require the successful occurence of the electroweak radiative 
symmetry breaking to restrict the parameter space. As a result the upper mass bound m^o 
is largely restricted compared with the one estimated without imposing this condition. 
We point out some features of m^o related to the initial value of h t and discuss why the 
models with more extra matters 5 + 5 of SU(5) could bring the larger maximum value of 
m h o. 
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The existence of a rather light CP-even neutral Higgs scalar is a common feature of the 
supersymmetric extension of the standard model |], ||. It is an important evidence of this 
kind of models and it is very useful to know its possible upper bound for the judgement 
of the consistency of the models. 

The next to the minimal supersymmetric standard model (NMSSM) is the simplest 
extension of the minimal supersymmetric standard model (MSSM) In this model a 
singlet chiral supermultiplet S is introduced and a fi term in the MSSM is replaced by 
a Yukawa coupling \SHiH 2 with the usual Higgs doublet chiral superfields Hi and H 2 . 
The superpotential of the Higgs sector in this model is expressed as 

H^NMSSM = XSHiH 2 + + " " " • (1) 

If the scalar component S of S gets a vacuum expectation value, the fi scale appears as 
X(S). In this model the \x problem in the MSSM is potentially solvable when the tree 
level \i term is forbidden due to a suitable symmetry j|. 

An interesting feature of this model is the fact that the mass bound of the lightest 
neutral Higgs scalar can be estimated with no dependence on the soft supersymmetry 
breaking parameters at the tree level f|, 

2X 2 
9i + 92 

where tan/3 = {H 2 )/{H^}. This bound is mainly controled by the value of tan/3 and 
the bound of A. Its dependence on the soft supersymmetry breaking parameters appears 
through the loop correction to the effective potential as a result of the large top Yukawa 
coupling h t 0. All of these effects are related through the renormalization group equations 
(RGEs). Many works have been done on this aspect || [/J. In Ref.0 it has been suggested 
that the additional extra matters such as 5 + 5 of SU(5) can heavily affect the mass bound 
by changing the running of gauge couplings, h t and A but not affecting the scale of the 
gauge coupling unification, which is a great success of the MSSM. 

In the NMSSM the radiative generation of nonzero (S) or the fi scale is a very im- 
portant aspect. From this point of view it seems to be rather crucial to study how the 
requirement of the successful occurence of the radiative symmetry breaking || affects 
on the upper mass bound of the lightest neutral Higgs scalar. In this note we investi- 
gate this problem by finding the radiatively induced minimum of the effective potential 
parameterized by (tan/3, (S)) in the models with extra matters n(5 + 5) of SU(5). 
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The one-loop effective potential due to the large top Yukawa coupling || |10 
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where Q is a renormalization point and m 2 is the eigenvalue of the stop mass matrix 
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The correction to Eq. (1) due to this one-loop effective potential is expressed by using 
these mass eigenvalues as 
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Here we used the potential minimum condition to eliminate the soft scalar masses of Higgs 
fields and took a field basis which is used to derive Eq. (2). The mass matrix in Eq. (4) 
depends on (tan/3, (S)) other than the Yukawa couplings and the soft SUSY breaking 
parameters. 

In order to fix this matrix we need to determine these values as the ones at the po- 
tential minimum. It is a nontrivial problem whether such values of (tan/3, (S)) can be 
radiatively realized starting from certain sets of the Yukawa couplings and the soft SUSY 
breaking parameters. Our task is to estimate m 2 h0 {= m^ 2 + Am^o) numerically for the 
parameter sets which can radiatively realize the phenomenologically acceptable potential 
minimum. We determine such parameter sets so as to satisfy the following conditions: 
(i) starting from the suitable initial values of parameters, the radiative symmetry breaking 
occurs successfully and the following phenomenologically required condition is satisfied at 
the potential minimum, 



(#i) 2 + (H 2 ) 2 = (174 GeV) 2 , m t = 174 GeV, (6) 

(ii) m 2 which corresponds to the one of diagonal elements of the 3x3 neutral Higgs mass 
matrix should be smaller than other two diagonal components ||, 

(iii) the experimental mass bounds on the charged Higgs bosons rriH±, charginos m x ±, 
stops fh tl and gluinos M3 are satisfied. These masses, except for M3, are dependent on A 
and (S). We require the following values for them: 

m H ± > 65 GeV, m x ± > 72 GeV m t2 > 67 GeV, M 3 > 173 GeV, (7) 
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(iv) the vacuum should be a color conserving one [11] . 



In this study we solve a set of RGEs which are composed of two-loop ones for di- 
mensionless couplings and one-loop ones for dimensional SUSY breaking parameters, for 
simplicity. As the initial conditions for the SUSY breaking parameters we take 

™% = (7*^) 2 , M a = M, A t = A K = A\ = A, (8) 

where m is the universal soft scalar mass. We introduce the nonuniversality represented 
by 7j only among soft scalar masses of Hi, Hi and S to make it easy to find the radiative 
symmetry breaking solutions. This will be taken as 0.8 < ji < 1.2. These initial condi- 
tions are assumed to be applied at the scale Mx where the coupling unification of SU(2)^ 
and U(l)y occurs. We donot require the regolous coupling unification of SU(3)c but only 
impose the realization of the low energy experimental value 0. The initial values of the 
parameters are surveyed through the following region, 

0<h t <1.2 (0.1), -2.0<k<0 (0.2), 0<A<3.0 (0.2), 

< M/M s < 0.8 (0.3), < m/M s , \A\/M S < 3.0 (0.5), (9) 

where in the parentheses we give the interval which we use in the survey of these pa- 
rameter regions.]] We also assume that the RGEs of the model are changed from the 
supersymmetric ones to the nonsupersymmetric ones at a supersymmetry breaking scale 
Ms, for which we take Ms = 1 TeV as a typical numerical value [H, [7|. 

To estimate the one-loop effect Eq. (5), it is necessary to know the values of (tan (3, (S)) 
at the potential minimum. If we impose the radiative symmetry breakling condition, such 
a potential minimum has to be realized as a result of RGEs solution. In order to see the 
effect of the radiative symmetry breaking condition on the upper mass bound m h o we 
calculate it under two situations. We impose the full condition of (i) in a case (I). On the 
other hand, in a case (II) we require only Eq. (6) but donot require that it is realized at 
the potential minimum^. We take the number of extra matters as n = 3. The perturbative 
unification of the gauge coupling requires n < 4. Later we will also discuss other cases. 

1 Since the sign of k and A affects the scalar potential, we need to investigate both sign of them. 
However, a negative n seems to cover almost solutions for the positive n. Here we give the only result in 

the case of the negative k. 

2 In other words, this case corresponds to the situation that Eq. (6) may be satisfied at the minimum 

which is obtained from the unnatural initial soft scalar masses. 
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Fig. 1 Scatter plots of the solutions for the case (I) and the case (II) in (tan/3, u(= (S)) and (tan/3, A) 
planes at the nit scale. In 1(a) only the solutions in (I) are drawn by circles. 

At first, we show how the radiative symmetry breaking condition restricts the param- 
eters strongly relevant to m^o. In Fig. 1 we give the plots of the solutions in (tan/3, (S)) 
and (tan /3, A) planes. Since we assume that only the top Yukawa coupling is large, the 
consistent tan/3 cannot be so large^J We take it as tan/3 < 15.0. In the case (II), (S) can 
take the value in such a wide range as 70 GeV < (S) < 45 TeV, which is not plotted in 
Fig. 1(a). These show that the value of A(rat) and (S) are heavily restricted in the case 
(I) compared with the case (II). We should remind that A(m t ) and (S) affect m^o through 
Eqs. (2) and (4). As seen from Fig. 1, the //(= \{S)) scale can take a value in a rather 
wide range to cause the successful radiative symmetry breaking. 

In Fig. 2 we show the allowed region of m h o in both cases for the corresponding 
parameter sets to the ones of Fig. 1. This shows that the imposition of the consistent 
occurence of the radiative symmetry breaking can strongly affect the estimation of m^o . 
The boundary value of m^o can be changed by a few percent to ten percent. We can 
also find some parameter dependences of m h o in these figures. In Fig. 2(b) we restrict 
the initial soft scalar mass as fh = 1 TeV. The larger values of soft scalar mass rh and 
(S) tend to realize the larger value of m^o. The value of (S) determines the off-diagonal 



3 To estimate tan/3 we take account of the translation of the pole mass to the running mass jlz 
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Fig. 2 Allowed regions of m h o. The boundary values are shown by the dashed and dash-dotted lines 
for ht > 0.5 and ht < 0.4, respectively. The displayed values of the top Yukawa coupling fit and the soft 
scalar mass fh are the ones at the unification scale Mx- 



component of the stop mass matrix (4). This shows that the larger stop mixing tends to 
make m h o larger. 

It is useful to comment on the existence of two branches of the solutions which show 
a very different behavior in Figs. 1 and 2. This is a common feature in both case (I) and 
(II) and also in the models with a different n. For simplicity, we take the case (II) as an 
example to discuss this feature here. We refer the one of the smaller tan/5 as a branch 
Bl and the one of the larger tan/? as a branch B2. They are divided by an initial value 
of ht as shown in Table 1. The interesting feature of two branches can be clearly seen in 
Fig. 1(b) and Table 1. The branch B2 corresponds to the smoothly extending solutions of 
X(m t ) to the large tan/3 region and the branch Bl comes from the ones which are confied 
in the small tan f3. The branch B2 has a similar lower bound of tan j3 and the similar 
maximum value of m^o ~ 143GeV at tan (3 > 9 for the different n values. If we note that 
at the larger tan (3 region the second term of Eq. (2) can be neglected and Eq. (2) reduces 
to the one of the MSSM, this behavior can be understood. On the other hand, the branch 
Bl strongly depends on n. There are solutions with a little bit larger maximum value of 
A at the smaller tan j3 according to the increase of n. This results in the larger maximum 
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Table 1 Extra matter effects and the branch structure in the case (II). 



value of m^o for the larger n. The reason can be mainly found in the (3 dependence of 
Eq. (2). Decreasing n, the maximum value of tan/5 of the branch Bl increases, where the 
maximum value of \{m t ) is realized. The branch Bl disappears finally at n = since it 
is difficult to realize 0.95 < h t (m t ) < 1.35 corresponding to 1 < tan/3 < 15 starting from 
the initial value of h t used here. We should remind that the value of h t (m t ) is strictly 
restricted by Eq. (6). In the case (I) we also find the similar qualitative feature discussed 
here in the case (II). 

Finally we want to present one comment. There is a rather big difference of the 
density of the radiative symmetry breaking solutions in the models with the different n. 
For example, in the present parameter setting the number of solutions in n = 4 is very 
smaller than the ones in n = 3. The finer tuning of marameters seems to be necessary to 
find the solution in the n = 4 case compared with the n = 3 case. 

In summary we studied the upper bound m h o of the lightest neutral Higgs scalar 
mass in the NMSSM using the RGE analysis. We required the successful occurence of 
the electroweak radiative symmetry breaking starting from the suitable initial values of 
parameters. This condition substantially constrains the allowable parameter space and as 
a result the mass bound m^o is heavily restricted compared with the one obtained without 
imposing this condition. We discussed the typical feature related to the initial value of h t 
and also why the larger n models could bring the larger maximum value of m h o. 

This work has been partly supported by the a Grant-in-Aid for Scientific Research 
from the Ministry of Education, Science and Culture(#11640267 and #11127206). 
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